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The partition function of the double-layer XY model in the (dual) Villain form is com- 



■ puted exactly in the limit of weak coupling between layers. Both layers are found to be 

cn , 

i locked together through the Berezinskii-Kosterlitz-Thouless transition, while they become 

O 

decoupled well inside the normal phase. These results are recovered in the general case 
! of a finite number of such layers. When re-interpreted in terms of the dual problems of 
lattice anyon superconductivity and of spin-liquids, they also indicate that the essential 
^ ! nature of the transition into the normal state found in two dimensions persists in the case 
' of a finite number of weakly coupled layers. 
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The layered crystal structure common to high-temperature superconductors is known 
to be reflected in their electronic properties.-^ There are suggestions in the theoretical lit- 
erature that the effect of strong repulsive interactions present within a single such copper- 
oxygen layer result in exotic anyonic superconducting states and/or spin-liquid states. 
Recently, the author has shown that both of these groundstates undergo a phase transi- 
tion into a normal phase at non-zero temperature, where the corresponding particle-hole 
excitations become deconfined upon heating through the transition.^ The transition itself 
is dual to that of the two-dimensional (2D) Coulomb gas.^'® In the present context of oxide 
superconductivity, therefore, one is naturally led to the question of how a finite number of 
such weakly coupled layers behave. 

In response to this query, we shall study the problem of a finite number of weakly 
coupled layers of square-lattice Xy-models in the Villain form,^ which is known to lead to 
a dual theory of (layered) compact quantum electrodynamics (QED) in the strong- coupling 
limit. The Abelian gauge-field that appears here describes chiral spin-fiuctuations in the 
case of the 2D spin-liquid state, while it represents the statistical flux-tube attached to 
each psuedo-fermion in the case of anyon superconductivity.^"^'-'^-'^ Korshunov has given an 
approximate solution of the XY-model with infinitely many weakly coupled layers, where 
he finds a Berezinskii-Kosterlitz-Thouless (BKT) transition within each layer, followed by 
an inter-layer decoupling transition at much higher temperature that is characterized by 
the binding of vortex rings (fluxons) lying in between consecutive layers. Note, therefore, 
that the lower-temperature superfiuid BKT transition involves the entire three-dimensional 
(3D) lattice in this scenario, which is consistent with previous Monte Carlo simulations. -"^^ 
In this Letter, we first solve the weakly coupled double-layer XY model in the Villain form 
exactly, recovering Korshunov's approximate results for the case of infinitely many layers 
in the processes. The intra-layer phase auto-correlation function of this double-layer XY- 
model is computed, where we find that it falls off algebraicly with distance below the BKT 
transition, as expected. Notably, the exponent of the correlation function in question 
vanishes linearly with temperature, yet the coefficient is only half as big as that of the 
single-layer case.^ This indicates that the existence of neighboring layers enhances long- 
range phase correlations, as expected. In addition, the inter-layer phase autocorrelation 
function is computed, where we find that it coincides with the prior intra-layer phase 
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autocorrelation function at temperatures below the decoupling transition. This, on the 

other hand, gives strong support to the claim that the double-layers become locked together 

at temperatures below the latter transition. Last, we show that the above double-layer 

results persist in the general case of a finite number of weakly coupled layers, suggesting 

that the superfluid transition exhibited by the XY model with a finite number of layers is 
3D-like.i2,i3 

When the above results are re-interpreted in terms of the dual problem of double- 
layer anyon superconductivity and spin-liquids, they yield that a confining string exists 
between particle-hole excitations at temperatures below the deconfinement transition. -"^^ 
Note that the consecutive layers have already become decoupled at this stage via the 
inter-layer fiuxon-antifiiixon unbinding transition that occurs at much lower temperature 
in this case. Taking into account arguments given previously by the author,^ this implies 
that the basic character of the single-layer particle-hole deconfinement transition into the 
normal state persists in the case of double-layer anyon superconductors and spin- liquids; 
i.e., the transition remains 2D. 

Let us first consider the generalized correlation function (exp[i = Z[p]/Z[0] 

of the A^-layered XY-model, where the partition function Z\p\ = J 'D(j)(r)e~^^^ has an 
energy functional 

N N-1 

T = E E E - cos[^M<^(^"'' 0]} +p±Yl E^i - ^^«['^(^"'' - m i+m 

1=1 f H=x,y 1=1 r 

N 

-iJ2Y.P(r,l)cj>ir,l). (1) 
1=1 f 

Here r ranges over the square-lattice and I denotes the layer index, while P\\^± = J\\^±/kBT 
is proportional to nearest- neighbor coupling-constants satisfying Jy 3> J±. Also, A^(/}{r) = 
4>{r + fi) — 0(r) is the lattice difference operator. The particular correlation to be probed is 
determined by the fixed integer field p(r, /); e.g., the choice p(±r/2, 0, 1) = ±1 corresponds 
to the intra-plane auto-correlation function. After making the usual (low-temperature) Vil- 
lain substitution for the exponential above,^ ^-Pli-cose] _^ (^2nP)-^/^ T.n=-oo e^^^e-^'/^/J^ 
and then integrating over the phase-field, we obtain the following dual representation equiv- 
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alent to AT-layered compact QED in the strong-coupling limit (modulo a const ant ):-'^° 

exp 



1 1 



1=1 f ^ 1=1 f 

(2 



where n^(r) is an integer link-field on the layered lattice structure of points r = (r, 
with fj, = x,y,z and n — {nx,ny). In the context of the dual systems of layered anyon 
superconductors and spin-liquids, n(r, I) represents the confining statistical electric field 
experienced between particle-hole excitations about such groundstates.^ Here one must 
make the identification I3^]<^ g\ ^^^fiuJa/kBT' , where y denote the corresponding gauge- 
field coupling constants, and where T' is the temperature in these dual systems. The Debye 
frequency scale, ujq = cqo"^, is set by the characteristic velocity cq, which represents the 
spin-wave velocity in the case of spin-liquids, and the zero-sound speed in the case of 
anyon superconductivity.^ Note that the case of weakly coupled planes, Jx <^ also 
corresponds to weakly electro-magnetically coupled layers in the dual system, where the 
coupling-constants satisfy <^ g'^'^. 

Double-layer. Consider now the dual representation outlined above for the special 
case of two coupled XY-model/compact QED layers. This problem may be solved ex- 
actly in the limit of weak interlayer coupling, J± <S Jy, by first making the decomposi- 
tion n(r, 1) = n'(r, 1) — n-{f) and n(r, 2) = n'(r, 2) -|- n-{f), such that the intra-layer 
fields n'{rj) satisfy V ■ n'\r = J2r' P(''^')^rr' , while the inter-layer field n-{f) satisfies 
V ■ n-\f — nzif), with V = (A^;, Ay). We now take the customary potential representa- 
tion n_ = — V$ for the inter-layer field, which yields $(r) = ^j^., G^'^\f—f')nz{r'), where 
G'(2)(^ = (27r)-2 d'^k{e^^-^-l)[4:-2cos{ka;a) - 2cos{kya)]-^ is the Greens function for 
the square lattice. Note that fluxon-charge neutrality is presumed;^ i.e., X^f '^zlr) = 0- 
The requirement that the original field n be integer valued indicates that the intra-layer 
field n' is not integer-valued, in general. However, in the presently considered limit of 
weak coupling, J±/J\\ — > 0, the concentration of inter-layer (fluxon) charge, is expo- 
nentially small. This implies that the inter-layer field, is equally small in magnitude, 
and that the intra-layer field, n', is very close to being integer valued. After making a 
suitable (lattice) integration by parts of the energy functional in Eq. (2), we thus obtain 
the factorization Z = Z^qZ^qZqg for the partition function in the limit of weak coupling 
between layers, where the intra-layer factors represent 2D discrete gaussian models that 



read 



{n'(r,l)} r' 



exp 



2/^11^ 



(3) 



and where the inter-layer factor is given by the 2D Coulomb gas ensemble 



ZcgIp]^ J2 expi--^^2K(r)-p(f,l)+p(f,2)]$(r)--i-^n2(r)i 

{n.(r)} I "^^11 r r J 

= I] expJ-^^2[n,(f)-Kr,l)+K^-',2)]G'(2)(f-fOn,(fO-2^X^ 

(4) 

Considering first the case of absent probe charges, p{r) = 0, we immediately see that each 
layer (3) undergoes a pure 2D Xy/discrete-gaussian model transition at ksTc ~ f Jy, 
while the inter-layer links nz{r) undergo an inverted 2D Coulomb gas binding transition 
at /cbT* = 47rJ|| in the limit of weak inter-layer coupling, J± <^ Jy. We therefore find 
that Korshunov's approximate results for the corresponding problem of infinitely many 
layers are exact in the present instance of double layers. Following this author, the latter 
high-temperature transition therefore corresponds to an inter-layer decoupling transition 
mediated by the binding of vortex rings lying in between the double-layers. 

Next consider the intra-layer phase autocorrelation function, Cii{R) = (exp[z0(0, 1) — 
1)]), for the XY-model, which corresponds to fixing two opposing unit charges 
separated by a distance R within the dual compact QED representation; e.g., we have 
p{ f, 1) = 6f^fi ~ R Pi^^ 2) = 0. By the previous factorization for the partion func- 
tion, we may therefore write Cii{R) = C\\{R)C±{R), where Ci\{R) = Z^l[p]/Z^l[0] and 
C±{R) = Zcg\p]/ZcgM- The "parallel" autocorrelation is given by that of the pure 
square-lattice XY model;^ i.e., C\\{R) = (ro/-R)''ii e~^/^ii , where the correlation length di- 
verges as /a ~ exp[^/(T/Tc — l)"*^/^] just above Tc,^ and where the exponent is given 
by ^11 = (27r/?||)~-^ at low-temperature and by ?7|| = 1/4 at Tg. Here, Inr^"^ = 1.6169. To 
compute the "perpendicular" factor C± (R) to the intra-layer correlation function, we first 
re-express the right side of Eq. (4) as Zcg[p]/ZcgM = exp{[G'[5j,(i?) - G^^\R)]/2i3^\}, 
where the meanfield renormalized Coulomb interaction is defined by^ 

{nz{r)} I " r,f' 



(r — r ') X 



X 



n 



(5) 



Therefore, for T < the above meanfield renormahzed Coulomb interaction is short- 
ranged, yielding C^(R) = exp[-G^^\R)/2p\\] ^ (i?/ro)^/^''^ii , while it is dielectrically 
renormalizedtoG'JJ^(i?) = e-^{T)G^^\R) at T > T*, yielding Ci(i?) ^ {R/ro)^^-'~'y^''(^\\ . 
Putting together all of these results, we obtain Cii{R) = (ro/-R)^e~^/^ii for the inter-layer 
correlator, where the temperature dependent exponent is given by 
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47r/3|| ' 



v= I; T^T, (6) 

Notice that the above correlations are strengthened by the second layer with respect to 
those of a single layer below the decoupling transition , while they approach those of an 
isolated single layer at high-temperature \\.im.T^oo^{T) = 1]. 

To compute the inter-layer phase auto-correlation function C\2{R) = (exp[z(/)(0, 1) — 
2)]), it becomes useful to consider the corresponding probe charges as a superposition 
of the previous intra-layer ones along with the minimum-distance inter-layer configuration 
p'(r, 1) = 5-^ and p'if^ 2) = r- The latter configuration can be incorporated into the 
inter-layer field n-{r), resulting in the modified constraint V ■ n-\r — '>^z{'f^ — ^. After 
appropriately modifying the manipulations beginning with Eq. (2), we ultimately obtain 
Ci2iR) = q|(i?)exp{-[G'(2)(i?) + G'^^(i?)]/2/?||} = (ro/i?)''i|-(i+^"')/^'^^iie-^/^ll for the 
inter-layer correlation function. Notice, therefore, that this inter-later correlator is equal 
to the intra-layer correlator, Cii(i?), at temperatures below the decoupling transition, 
T < T:^, where e~^{T) = 0. We see explicitly, therefore, that the double-layers become 
locked at temperatures below T*. 

In terms of double-layer compact QED, which describes the gauge-field dynamics in 
both double-layer anyon superconductors and spin- liquids,^ the above discussion implies 
that there exists a deconfinement transition at A;r hiOQ^ below which both intra- 

layer [Cii(i?)] and inter-layer [Ci2(-R)] particle-hole excitation in either system are confined 
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{n{k, 1) • n(-^,2)) = • n_(-fc)) = - „ , , , (7) 



by a statistical electric flux tube^° with string tension 7(T') = ksT' /^\\{T'). Note that 
this result is identical to that of a single layer. The present double-layer case, however, 
also shows a decoupling transition at ksTl = {An)~^g'^^huJo ~ /csT^'/lO between the corre- 
sponding in-plane statistical electric fields, n{f, I). In particular, the Fourier transform of 
the inter-layer correlation function for this field is given by 

A' 

where ^_l(T') denotes the Debye screening length for the 2D Coulomb gas of fiuxons. 
Above, we have used the potential representation ?2_ = — V$ for the "difference" field, 
along with the identity ($(^)$(— ^)) = ^(/c^ -|- for the Coulomb gas correlation 

function. Hence, the homogeneous inter-layer correlation (n(0, 1) • n(0,2)) jumps from 
—eQ^{Tl) = — |- to zero upon heating through decoupling transition-temperature T^. This 
indicates that the deconfinement transition at > Tl is 2D, since the layers have already 
become decoupled by then. Finally, we note that the specific-heat of the present double- 
layer model for strong-coupling compact QED naturally decouples into C„ = 2Cdg + C'cG7 
where Cdg and Ceo denote the contributions due to an isolated layer (3) described by 
the discrete-gaussian model^ and those due to the inter-layer Coulomb-gas ensemble (4), 
respectively. Given that Ccg(2^') shows a smooth peak just above^ T^, while Cdg(7") 
exhibits one just below^ ~ lOT^, we expect that the total specific-heat of the double- 
layer, Cv{T'), has two smooth peaks in the regime Tl <T' < T^. 

N -layers. Let us now generalize the previous discussion to the case of N weakly 
coupled XY models (1). To simplify matters, let us also consider the partition function 
(2) in the absence of probe charges, i.e., p = 0. The dual in-plane field can then be 
expressed as n(r, /) = n'(f, /) — n-{f^ I) +n-{f, ^ — 1), along with constraints V-n' = and 
V • n_ = Hz- Again, taking the customary potential form for the "difference" field, n_ = 

— * 

— V$, and making suitable (lattice) integration by parts, one obtains the factorization 
Z = ZcG^iLi^ for the partion function in the limit J±/J\\ — > 0, where the Coulomb 
gas portion (4) now has the generalized form 

r TV iV-l 

^cG= E exp -— EE[^-(^"''^-i)-^-(^"''0]'-^EE^'(^"''0 

{n^{f,l)} I ^'1 1=1 r 1=1 r 

f 1 

= E e^Pl E Et^-^^"' ^ - 1) - l)]G^^Hr-r')[n.{f\ I - I) - n^if', I)] 

{n;,{f,l)} I 1=1 r,f' 



with the fields at the boundary layers set to nz{r, 0) = = n^(r, A^). The key difference 
between the present Coulomb gas ensemble and the elementary one associated with the 
double-layer case (4) is that consecutive vertical fields, ^^(r, / — 1) and nz{f,l), of equal 
value make no contribution to the energy functional. However, in the limit of weakly cou- 
pled planes, Eq. (8) indicates that the concentration of dipolar (fiuxon) charge excitations 
between consecutive planes is exponentially low. This observation suggests that each set 
of consecutive layers decouple precisely like the isolated double-layer at T*, in agreement 
with the results obtained by Korshunov in his approximate treatment of the corresponding 
problem of infinitely many layers. At temperatures below T*, we expect that consecutive 
layers will be locked, with each layer in turn attaining quasi long-range order at tempera- 
tures below Tc ~ T^/10. In this picture, therefore, the superfiuid BKT transition involves 
the entire layered structure, making it 3D-like.^^ 

The above picture may fail, however, in the regime of moderately coupled to strongly 
coupled planes, where the assumption that dipolar (fiuxon) charge excitations be dilute is 
no longer valid. This is certainly true for the case of the isotropic problem with infinitely 
many layers (3D XY model), where the decoupling transition necessarily must be driven 
to infinite temperature due to the occurrence of a single 3D phase transition in the system. 
Hikami and Tsuneto have argued that vortex rings traversing many layers have a mini- 
mum mean radius given by r ~ iJ\\/ J±y^'^ in units of the lattice constants. Therefore, 
comparison of the latter to the total number of layers, A^, yields a minimum anisotropy 
condition J\\/J± ~ N"^ that insures the existence of the BKT superfiuid transition char- 
acteristic of weakly coupled double-layers. This indicates, however, that the XY model 
with infinitely many layers may not be in the same universality class as that of the weakly 
coupled double-layer, contrary to claims in the literature. 

In conclusion, we find that the superfiuid transition exhibited by a finite number (two 
or greater) of weakly coupled square-lattice XY model layers in the dual Villain form is 
3D-like^^ while remaining in the same universality class as that of the 2D XY model. 
This result suggests that the latter transition is driven not by the unbinding of isolated 
vortex-antivortex pairs in a single layer, which are in fact (linearly) confined, ■'^^'■'^^ but by 
the unbinding of lines of vortex-antivortex pairs passing perpendicularly through all of the 
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layers. We also find an inter-layer decoupling transition in the normal state at much higher 
temperature.^^ Last, in terms of the dual problems of lattice anyon superconductors and 
spin-liquids,^ these results indicate that the essential nature of the respective supercon- 
ductor/normal and spin-liquid/paramagnetic transitions found in 2D is preserved in the 
present case of weakly coupled layered structures. 
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Science Foundation grant DMR-9322427. 



9 



References 



1. The Physical Properties of High- Temperature Superconductors, vol. 2, edited by D.M. 
Ginsberg (World Scientific, Singapore, 1990). 

2. P.W. Anderson, Science 235, 1196 (1987). 

3. R.B. Laughlin, Phys. Rev. Lett. 60, 2677 (1988); R.B. Laughlin, Science 242, 525 
(1988). 

4. J.P. Rodriguez and B. DouQot, Phys. Rev. B45, 971 (1992). 

5. N. Read and S. Sachdev, Phys. Rev. Lett. 62, 1694 (1989); J.P. Rodriguez, Phys. 
Rev. B 41, 7326 (1990); A. Sokol and D. Pines, Phys. Rev. Lett. 71, 2813 (1993). 

6. J.P. Rodriguez, Los Alamos Report # LA-UR-94-3002 (1994). 

7. S.T. Chui and J.D. Weeks, Phys. Rev. B 14, 4978 (1976). 

8. P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987). 

9. J.V. Jose, L.P. Kadanoff, S. Kirkpatrick and D.R. Nelson, Phys. Rev. B 16, 1217 
(1977); C. Itzykson and J. Drouffe, Statistical field theory, vol. 1, chap. 4 (Cambridge 
Univ. Press, Cambridge, 1991). 

10. A. Polyakov, Phys. Lett. 72B, 477 (1978); L. Susskind, Phys. Rev. D 20, 2610 
(1979). 

11. X.G. Wen, F. Wilczek, and A. Zee, Phys. Rev. B 39, 11413 (1989). 

12. S.E. Korshunov, Europhys. Lett. 11, 757 (1990). 

13. S.T. Chui and M.R. Giri, Phys. Lett. A 128, 49 (1988). 

14. S. Hikami and T. Tsuneto, Prog. Theor. Phys. 63, 387 (1980). 

15. L.L Glazmann and A.E. Koshelev, Zh. Eksp. Teor. Fiz. 97, 1371 (1990) [Sov. Phys. 
JETP 70, 774 (1990)]; V. CataudeUa and P. Minnhagen, Physica C 166, 442 (1990); 
P. Minnhagen and P. Olsson, Phys. Rev. B 44, 4503 (1991). 



10 



